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The application of the Kudryashov method for finding exact solutions of the high order non- 
linear evolution equations is considered. Some classes of solitary wave solutions for the families of 
nonlinear evolution equations of fifth, sixth and seventh order are obtained. The efficiency of the 
Kudryashov method for finding exact solutions of the high order nonlinear evolution equations is 
demonstrated. 
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G\ ■ 1 Introduction 

The powerful and effective method for finding exact solutions of nonlinear ordinary differential 
equations was proposed in work [T]. In works author applied this method to construct the 

exact solutions of the nonlinear nonintegrable equations. The first modification of this method was 
presented in work [5 . The final and the most successful modification of this method was proposed 
in [6]. Thus, we refer this method as the Kudryashov method. The Kudryashov method allow us in 
the straightforward manner to construct solitary wave solutions for a wide class of nonlinear ordinary 
differential equations. The main idea of the Kudryashov method is to use special form of the singularity 
manifold in the truncation method [IH5] . This approach allows us to reduce the problem of constructing 
exact solutions to solving the overdetermined system of algebraic equations. The main advantage of the 
Kudryashov method is that we can more effectively construct exact solutions of high order nonlinear 
evolution equations in comparison with other methods for finding exact solutions [71113) . 

Using this method in works [IH3] exact solutions of the generalized Kuramoto-Sivashinsky equa- 
tion, the Burgers-Korteweg-de Vries eqution, the Bretherton eqaution and the Kawahara equation 
were obtained. Traveling wave solutions for class of third order nonlinear evolution equations were 
constructed in [14] with help of the Kudryashov method. The Kudryashov method was used for 
constructing traveling wave solutions of several nonlinear evolution equation in works [15H16] as well. 

The aim of this work is to demonstrate efficiency of the Kudryashov method for finding exact 
solitary wave solutions of high order nonlinear evolution equations. For this purpose we consider three 
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families of nonlinear evolution equations of fifth, six and seven order. This families of equations have 
some applications in physics and other fields of science. We present classes of solitary wave solutions 
for considered families of nonlinear evolution equations. 

This work is organized as follows. In the next section we give brief description of the Kudryashov 
method algorithm. In the sections 3,4 and 5 we construct solitary wave solutions for the families of 
nonlinear evolution equations of fifth, six and seven order respectively. In the last section we summarize 
and discuss our results. 



2 Method applied 

The aim of this section is to present the algorithm of the Kudryashov method for finding exact 
solutions of the nonlinear evolution equations. To reach this purpose we will follow the works [5U5UT3]. 
Let us consider the nonlinear partial differential equation in the form 

Ei[u u u x , ...,x,t] = (1) 

Using the following ansatz 

u(x, t) — y(z), z — kx ~ ujt. (2) 
from Eq. (fTJ) we obtain the ordinary nonlinear differential equation 

E 2 [-ujy z ,ky z ,k 2 y zz ,k 3 y zzz , . . .} = 0. (3) 

Now we show how one could obtain the exact solution of the Eq. ([3]) using the approach by 
Kudryashov. This method consist of the following steps [SJE1H3] • 
The first step. Determination of the dominant terms. 
To find dominant terms we substitute 

y = z\ (4) 

into all terms of Eq. (j3|) . Then we compare degrees of all terms in Eq. ([3]) and choose two or more with 
the smallest degree. The minimum value of p define the pole of Eq. ([3]) solution and we denote it as 
N. We have to point out that method can be applied when N is integer. If the value TV is noninteger 
one can transform the equation studied and repeat the procedure. 

The second step. The solution structure. 

We look for exact solution of Eq. ([3]) in the form 

y = a a + axQ(z) + a 2 (z)Q{z) 2 + ... + a N Q{z) N , (5) 
where et^ - unknown constants, Q(z) is the following function 

^ = ih- (6) 

This function satisfies to the first order ordinary differential equation 

Q z = Q 2 - Q. (7) 

The Eq. ([7]) is necessary to calculate the derivatives of function y(z). 
The third step. Derivatives calculation. 
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We should calculate all derivatives of function y. One can do it using the computer algebra 
systems Maple or Mathematica. As an example we consider following case: 

The derivatives of junction y{z) in the case of N = 2 can be written in the form 

y = a + aiQ + a 2 Q 2 , 
y z = -a 1 Q+(a 1 -2a 2 )Q 2 + 2a 2 Q 3 , (8) 
Vzz = aiQ + (4a 2 - 3 ai )Q 2 + (2oi - 10a 2 )Q 3 + 6a 2 Q 4 . 

The relations JHJ) can be generalized for any value of N. Differentiating the expression ([5]) with 
respect to z and taking into account we have 

N 

Vz = ^<M(Q - 1)Q\ 

1=1 (9) 

i=l 

The high order derivatives of function y(z) can be found in works 5,6 . 
The fourth step. Defining the values of unknown parameters. 

We substitute expressions © in Eq. ((T|). After it we have take into account (JSJl. Thus Eq. fl} 
takes the form 

P[Q{z)] = 0, (10) 

where P[Q{z)\ - is a polynomial of function Q(z). Then we collect all items with the same powers of 
function Q{z) and equate this expressions equal to zero. As a result we obtain algebraic system of 
equations. Solving this system we get the values of unknown parameters. 

The Kudryashov method is a very powerful method for finding exact solutions of the nonlinear 
differential equations. It has a set of advantages. They are: 

1. The first step is not necessary, because all redundant terms in ([5]) becomes equal to zero when 
we start to define unknown parameters from the system of algebraic equations; 

2. We construct the solution as a set of Q-blocks. The function Q does not contain any parameters. 
It is very comfortable, because all of them is in equation; 

3. Given method can be easily programmed in Maple or Mathematica, because we use the substi- 
tutions and (JHJ) in Eq. © it takes the polynomial form; 

4. This method is powerful and effective even if we construct the exact solutions of the high order 
nonlinear evolution equations. 

5. It easy to show that tanh, coth, (G'/G) methods and Kudryashov method can be reduced to each 
other [B]. Moreover the Kudryashov method gives the same results as an Exp- function method. 
However it is well known that Exp-function method can not be applied for the equations of hight 
order; 

Let us give several examples to demonstrate it's efficiency. 
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3 Exact solitary wave solutions of the fifth order evolution 
equation 

As an example let us consider the fifth order nonlinear evolution equation in the form 

u t + uu x + lQuu xxx + 20u U x ~\~ OLli X x ~\~ PUxxx~\~ 

(11) 

~\~ r )'U>xxxx ~\~ Uxxxxx — 0. 

Eq. (jlip is new and does not present in the periodic literature. However it is an interesting equation 
because this equation consist of the fifth order Korteweg-de Vries equation with additional dispersive 
and dissipative terms. The equation (jf f p can arise in the physical applications when we consider the 
wave processes in active dispersive-dissipative media. 

Using the traveling waves ([2]) we have 



-wy z + kyy z + 10k yy zzz + 20k y z y zz + 30ky y z + ak y zz + f3k y zzz + 

y zzzz + k y ZZ zzz — 0. 



(12) 

The pole of the Eq. (| 12[) is equal to N = 2, thus we look for exact solution in the form 

y = a + aiQ + a 2 Q 2 , (13) 

where ao, a\ and a 2 - are unknown constants. 

Substituting © in Eq. (fT2")l and taking into account (| 13[) we obtain the polynomial of function 
Q(z). Collecting all terms with the same power of function Q(z) and equate this expressions to zero 
we obtain the system of algebraic equations. Solving this system we find that solution of Eq. (IT2"j) 
exists only in eight cases. They are 



fc (1)2) =±-^ fc (3 , 4 ) = ±^, fc(5, 6) =±-V988 7 2 + 2646/3 -441, 

21 7 84 (14) 

fc (7)8) = ± — V5 7 (972 7 3 + 343a). 

However in the case of k = fc( 5 6 ) and k = ktr,8) t ne parameters and solution presentation is very 
cumbersome and we do not give them in the present manuscript. 

In the case of k = fc(i,2) an d k = fc(3,4) the values of parameters a, /3, w, ao, a±, a 2 are 

7 5 5 44 3 

fc =21' a= 84 7 - 14 7/? -441 7 ' 
P , 16 7 2 1 4 2 2 2 

fl0 = 20 + 2205 ~40' ai= 147 7 ' ° 2 = " 147 7 ' (15) 
7 (9472 7 4 - 10416 7 2 - 9261 - 37044/3 + 62496/3 7 2 + 111132/3 2 ) 



31116960 

7 5 5 44 3 

k= -YV a= 84 7 - 14 7/3 -441 7 ' 
(3 46 7 2 1 2 2 

an = 1 . fll =0, 0.2 = 7 , 16) 

U 20 2205 40' '147 ' ' V ; 

7 (9472 7 4 - 10416 7 2 - 9261 - 37044^ + 62496/3 7 2 + 111132/3 2 ) 

31116960 ' 
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ai = 0, a 2 
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1507 60 1225 49 
67228a 2 - 7056a 7 3 + 347328 7 6 - 85757 2 
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24010007 



a 



7a 



1 



7 ' 

228 



:7 j 01 



1 18 

6 ~ 7 
36 



'7 



:7 , a 2 



18 



1507 60 1225 ' 49 ' 49 

67228a 2 - 7056 7 3 a + 347328 7 6 - 85757 2 
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24010007 ' 
The solutions of the Eq. (fTS"]) which corresponds to the relations (|T5|) — (fT8|) are 



y(z) = a + — 1 2 [2-Q(z)]Q(z), 
2 7 2 

y(z) = a - —Q{z) 2 , 

y{z) = 00 - ^g-Oi 2 ) » 
y (z) = a + Hi! [2 - Q(z)] Q(z). 

Graphical presentation of solutions (|T9")) . (|2T)|) is shown on Fig. [TJ 




0.041 
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Figure 1: Exact solutions (fT9| and ([20| - 1, 2 respectively at (3 = 7 = 1. 



(17) 



(18) 



(19) 
(20) 
(21) 
(22) 



The solutions (|21j) . (|22j) have the same structure so we have decided not to picture them. 



4 Exact solitary waves of the sixth order evolution equation 

Let us consider the following equation 



u t + uu x + au xx + 0U 



= 0. 



(23) 



This equation was proposed in work |17j for describing the longitudinal seismic waves in a vis- 
coelastic medium. Moreover this equation is used for modeling the "soft" - type of the turbulence [T5] . 
Also this equation describes the chemical reactions in reaction-diffusion systems [19] . The numerical 
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modeling of the wave processes describing by Eq. ([23]) was performed in work [3D] . Exact solutions of 
this equation was obtained in work |21j . 

Taking into account the traveling waves <j2j) in Eq. (|23|) we obtain 

- ujy z + kyy z + ak 2 y zz + /3k 4 y zzzz + jk 6 y zzzzzz = 0. (24) 

Integrating the Eq. (j2~4")l we have 

y 2 

Ci-iuy + k— + ak 2 y z + Bk 4 y zzz + -/k 6 y zzzzz = 0. (25) 

Dominant terms of the Eq. (f2"5"]) are k 6 y zzzzz , ky 2 /2. Thus, the pole order of the Eq. (|2~5"|) solution 
is N = 5. So we look for solution in the form 

y(z) =a a + aiQ + a 2 Q 2 + a 3 Q 3 + a 4 Q 4 + a 5 Q 5 , (26) 

where ao - as - constants to be determined. 

Using the in Eq. (|2"5"j) and taking into account ansatz we obtain a system of algebraic 
equations. Solving this system we hnd four real families of unknown parameters. However we give 
only three because one of them is cumbersome. The following families are 



114 ~ 8 lo 2 127008fc 7 ^ 2 
— k ^ 7 = "55F' ° 1= 2k 121— 



uj 504fc 3 /3 „ 10080fc 3 /3 

ao = -, 7-: — , ai=a 2 = 0, a 3 = — , (M) 

k 11 11 

15120fc 3 /3 6048fc 3 /3 
a 4 = — , «5 = ^— 

219, 2 a 5/3 w 2 99225& 7 /3 2 

a =^ k ^ ^=~u¥> Cl= 2k~ 968 ' 

u 315/c 3 ^ 19845fc 3 /3 „ , oc ^ 

fl o = T + ^T^> a i= ' «2 = 7i — ' «3 - 6930£; 3 /3, (28) 

K 22 11 

94500/c 3 / 3 37800fc 3 ^ 
a 4 = — , «5 = ^— 

3259 „. 2 1 /3 _ w 2 321489fc 7 /3 2 

a = _ W*' 7 = ~iW Cl= 2k 242^' 

u 567fc 3 /3 „ 1890fc 3 /3 3780fc 3 /3 fom 

a o = j — — ' ai = ' 0,2 = — n — ' a3 = — n — ' 

7560fc 3 /3 3024fc 3 /3 
a 4 = jj— , "5 = ^— 

From relations (|27|) . (|28|) and ([29)) we see that the real solution of Eq. (|25|) exists in the case of 

2 2 = 275 

114/3' P 114 7 ' 

,„ 44 a „ 9 1095 , s 

k = , B 2 = cry, 30) 

219/3' M 1936 h { ' 

2 _ HO a 2 _ 3259 

~~ 3259 0' ^ _ T2T00 a7 - 
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The solutions of Eq. (f2"5"j) which corresponds to the relations (|27p , (f!?5|) and takes the form 
uj Bk 3 

y(z) = - + — (-504 + 10080Q(z) 3 - 15120Q(z) 4 + 6048Q(z) 5 ) , (31) 
k 11 ' 



y(z) = ~ + — ( ^ - 19845Q(z) 2 + 76230Q(z) 3 - 94500Q(z) 4 - 
k 11 y 2 

+37800Q(z) 5 ) , 

uj Bk 3 

y(z) = - + — (-567 + 1890Q(z) 2 + 3780Q(z) 3 - 7560Q(^) 4 + 
k 11 

+3024Q(z) 5 ) . 



(32) 



(33) 



The solutions (|3ip. (|32p of the Eq. (|2"3")l is presented on Fig. [5J The third solution has the form of 
the kink as well as (|3ip. 
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Figure 2: Exact solutions (|3~Tj) and (|32| - 1, 2 respectively at B = uj = k = \. 



5 Exact solitary waves of the seventh order evolution equation 



Let us consider the nonlinear evolution equation of seventh order 

~l~ ^ ^a; ~l~ OlU xxx -\- fiUxxxxx ~t~ ^^xxxxxxx 0, 



(34) 



where n = 1, 2, 3. 

For the first time this equation was considered in work [5] in the case n = 2. We generalizes the 
results of this work. 

Taking the traveling wave ansatz u{x, t) = y(z), z = kx — wt into account from Eq. (l34p we have 



- uy z + ky n y z + ak 3 y zzz + /3k 5 y zzzzz + -fk 7 y. 
Integrating the Eq. (13"51 with respect to variable z we obtain 



zzzzzzz 



0. 



,n+l 



C\-ujy + k — — + ak 3 y zz + /3k 5 y zzzz + -fk 7 y z 
n + 1 



(35) 



(36) 
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In the case of n = 1 we have to look exact solution of Eq. (|3"6")l in the form 

y(z) = a + a\Q + a 2 Q 2 + a 3 Q 3 + a^Q 4 + a 5 Q 5 + a§Q 6 , 

because the pole order of Eq. (155)) solution is N = 6. 

Substituting (J9j) in Eq. (l36t with n = 1 and taking (|37|) into account we found the 
solutions 

y(z) = %■ + 7 fc 6 (-7100 + 166320Q(z) 2 + 332640Q(z) 3 - 
-1829520Q(z) 4 + 1995840Q(z) 5 - 665280Q(z) 6 ) , 

2/0) = r + 7 fc6 (-3600 + 665280Q(z) 3 - 1995840Q(z) 4 
fc 

+1995840Q(z) 5 - 665280Q(z) 6 ) , 
The corresponding parameters of Eq. and (f3"7) are 



(37) 

following 
(38) 

(39) 



/3 = -100fc 2 7 , a = 2159fc 4 7 , C x = — - 25205000fc 13 7 2 , 

2k 

a = ^ - 7100fc 6 7 , ai =0, a 2 = 166320fc 6 7, a 3 = 332640Ac 6 7 , 
a 4 = -1829520fc 6 7, a 5 = 1995840fc 6 7 , a 6 = -665280/c 6 7, 



(40) 



/3 = -50fc 2 7 , a = 769/c 4 7, d = — - 6480000fc 13 7 2 , 

2k 

a = % - 3600fcS, at = a 2 = 0, a 3 = 665280fc 6 7 , 



a 4 = -1995840/c 6 7, a 5 = -a 4 , 
where /3/ 7 < 0. Thus, the solutions (|3"8")l and (|39l) exist when 

,2 _ ^2 _ 1000 



100 7 ' r 2159 
. a _ /3 fl2 _ 2500 



a6 = -03, 



07, 



5O7' 



769 



07. 



The illustration of ([3"5f is given on Fig. [3] 
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Figure 3: Exact solution j38]) at w = 10, 7 = 0.001, fc = 1. 

The solution (l39l) has the same form. 



(41) 



(42) 



Let us consider the case of n = 2. The pole order of Eq. (j3"6")l solution in that case is equal to 
N = 3. Thus, we look for exact solution of Eq. (j3"6")l in the form 



= a + aiQ + a 2 Q 2 + a 3 Q 3 



(43) 



where ao, ai, a,2, 03 the constants to be determined. 
Taking into account ^ and relations (l4"51 in Eq. 



we have the solution in the form 



y(z) = ±6v/^105^fc 3 (1 - 6Q(z) 2 + AQ(zf) , 



where 7 < 0. 

Values of (3, a, w, C\ are defined by relations 



P = -83/c 2 7, a = 946fc 4 7 , w = -1260jk 7 , d = 0, 
a = ±6v/-1057fc 3 , ai =0, a 2 = ^36^-1057^, 



a 3 = ±24 v /-105 7 fc 3 , 
So, the solution (|4^|) exist when the wave number satisfies to equivalence 



83 7 ' P 946 T 



The solution flU) is a kink, see Fig. gj 

80 
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Figure 4: Exact solution (|44|) at 7 = — 1, A; = 1. 



(44) 



(45) 



(46) 



In the case of n = 3 the pole order of the Eq. (|36j) solution is = 2. So, we use the following 
ansatz 

y(z) = a + aiQ + a 2 Q 2 , (47) 

where ao, 01,02 - unknown constants. 

Then we substitute ([9]) and (|47|) in Eq. (5.3) with n = 3. Collecting all terms with the same power 
of function Q{z) and equate them to zero we obtain the system of algebraic equations on unknown 
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parameters. Solving this system we have 



5/3 2 

a = 14fc 4 7 H , w 



Ci 



ao 



28 7 
^315 k 
65856 7 8 / 3 



11767' 



(3920&V - 1764/3 7 2 fc 4 - 5^ 3 ) , 



(5^ 4 - 137200fc 8 7 4 + 9408/c 6 7 3 /3 + 392fc 4 7 2 /3 2 ) , 



-42-^73 3 ' a i= 4 V3157fc, a 2 



-ax. 



The solution of Eq. (j3"6")l with n = 3 which corresponds to (|4T>1) takes the form 

»(*) = -^^3 " (1 ~ WOW + • 

This (|49|) has the form of the traveling wave (Fig. [5]). 




(48) 



(49) 
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Figure 5: Exact solution (|49|) at/3 = 7 = A; = l. 



6 Conclusion 

In this work we have demonstrated efficiency of the Kudryashov method for finding exact solutions 
of high order nonlinear evolution equations. We have obtained solitary wave solutions for the three 
families of nonlinear evolution equations of fifth, six and seven orders. Graphical representation of 
these exact solution is presented. We believe that some of exact solutions obtained in this work are 
new. 

7 Acknowledgements 

This work was supported by the federal target programm "Research and scientific-pedagogical 
personnel of innovation in Russia" on 2009-2011. 



References 

[1] Kudryashov N.A. Exact soliton solutions of the generalized evolution equation of wave dynemics. 
J. Appl. Math. Mcch. 52 (1988) 361-365. 



10 



[2] Kudryashov N.A. Exact solutions of the generalized Kuramoto-Sivashinsky equation. Phys. Lett. 
A. 147 (1990) 287-291. 

[3] Kudryashov N.A. On types of nonlinear nonintegrable equations with exact solutions Phys. Lett. 
A. 155 (1991) 269-275. 

[4] Kudryashov N.A. Singular manifold equations and exact solutions for some nonlinear partial 
differential equations. Phys. Lett. A. 182 (1993) 356-362. 

[5] Kudryashov N.A. Analitical theory of nonlinear differential equations. Moskow - Igevsk. Institute 
of computer investigations: 2004 [in Russian] 

[6] Kudryashov N.A. On one of methods for finding exact solutions of nonlinear differential equations. 
larXiv:1108.3288l l [nlin.SI] 16 Aug 2011 

[7] Parkes E.J. Exact solutions to the two-dimensional Korteweg-de Vries-Burgers equation. J. Phys. 
A: Math. Gen. 27 (1994) 27 L497-L501. 

[8] Parkes E.J., Duffy B.R. An automated tanh-function method for finding solitary wave solutions 
to non-linear evolution equations. Comput. Phys. Commun. 98 (1996) 288-300. 

[9] Malfliet W., Hereman W. The tanh method: I. Exact solutions of nonlinear evolution and wave 
equations. Phys. Scr. 54 (1996) 563-568. 

[10] Fan E. Extended tanh-function method and its applications to nonlinear equations. Phys. Lett. 
A. 277 (2000) 212-218. 

[11] Biswas A. Solitary wave solution for the generalized Kawahara equation. Appl. Math. Lett. 22 
(2009) 208-210. 

[12] Vitanov N.K., Dimitrova Z.I., Kantz H. Modified method of simplest equation and its application 
to nonlinear PDEs. Applied Mathematics and Computation 216 (2010) 2587-2595 

[13] Kudryashov N.A., Loginova N.B. Extended simplest equation method for nonlinear differential 
equations. Applied Mathematics and Computation 205 (2008) 396-402 

[14] Ryabov P.N. Exact solutions of the Kudryashov-Sinelshchikov equation. Applied Mathematics 
and Computation 217 (2010) 3585-3590. 

[15] Kabir M. M., Khajeh A., Abdi Aghdam E., Yousefi Koma A. Modified Kudryashov method for 
finding exact solitary wave solutions of higher-order nonlinear equations. Mathematical Methods 
in the Applied Science, doi: 10.1002/mma.l349 

[16] Gepreel K.A., Omran S., Elagan S.K. The Traveling Wave Solutions for Some Nonlinear PDEs in 
Mathematical Physics. Applied Mathematics 2 (2011) 343-347. 

[17] Nikolaevskii V.N. Dynamics of Viscoelastic Media with Internal Oscillators, in Recent Advances 
in Engineering Sciences, Ed. by S.L. Koh and C.G. Speciale (Springer, Berlin, 1989), 210-221. 

[18] Tribelsky M.I., Tsuboi K. New Scenario for Transition to Turbulence? Phys. Rev. Lett. 76 (1996) 
1631-1634 . 

[19] Tanaka D. Chemical Turbulence Equivalent to Nikolaevskii Turbulence, Phys. Rev. E 70 (2004) 
015202-1-015202-4 



11 



[20] Kudryashov N.A., Migita A.V. Periodic structures developing with account for dispersion in a 
turbulence model. Fluid Dyn. 42 (2007) 463-471 

[21] Kudryashov N. A., Nonlinear differential equations with exact solutions expressed via the Weier- 
strass function. Zeitschrift fur Naturforschung 59 (2004) 443-454. 



